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I. INTRODUCTION 


Electromagnetic boundary-value problems can be solved exactly via 
classical separation-of-variable analysis only for a few geometries, 
such as the sphere, spheroid, circular cylinder, elliptical cylinder, 
strip and wedge. Such solutions can be expressed in terms of a summation 
of a set of eigenfunctions which can be evaluated with a high speed 
computer. However, for problems involving complicated geometries the 
exact solution is not available. Thus, the approximate and numerical 
methods are of great significance. The research discussed herein is 
directed toward the numerical solution of the problems of radiation and 
scattering from conducting bodies of arbitrary shape. 

Two methods are available for electromagnetic modeling of continuous 
conducting surfaces with arbitrary shape: the. wi re-gri d model [1] and 

the surface-current model [2,3] using rectangular-pulse bases. Both 
methods have similar limitations with the maximum cell width restricted 
to approximately A/10. Unless the conducting body is symmetric or is 
a figure of revolution, computer storage requirements have limited the 
moment-method application to bodies with surface area not exceeding one 
or two square wavelengths. 

In this report, the piecewi se-sinusoi dal reaction technique applied 
by Richmond [4,5] to thin-wire antennas is extended to scattering by 
conducting bodies of arbitrary shape. The continuous conducting surface 
is divided into cells and the surface-current distribution is expanded 
in overlapping sinusoidal bases. Via an application of Galerkin's 
method, the integral equation formulated wi th the zero-reaction concept [6] 
is reduced to a matrix equation. From a physical viewpoint, this 
variational solution follows from enforcement of reaction tests with 
an array of sinusoidal electric test sources. The current distribution 
over the conducting surface is determined via matrix inversion. 

Finally the scattered field is obtained by integrating the surface 
currents. 

The remaining text presents the detailed theory of the sinusoidal 
reaction formulation for radiation and scattering problems. The 
time dependence eJ^t f 0 r the time-harmonic source is understood and 
suppressed. Numerical results are presented for the scattering patterns 
of rectangular plates and Corner reflectors and the radiation patterns 
of corner-reflector antennas. 

The experimental data on corner-reflector antennas were measured by 
Melvin C. Gilreath at NASA Langley Research Center. We appreciate 
sincerely his kind permission to reproduce these antenna patterns. 

II. THE REACTION TECHNIQUE 

The reaction concept and its applications have been discussed 
by Rumsey[6], Cohen[7], Harrington[8] and Richmond[9]. 


Consider the exterior scattering problem illustrated in Fig. la. 
In the presence of a dielectric or conducting body, the impressed 
electric and magnetic currents (Jj ,Mj ) generate the electric and 
magnetic field intensities (jE,H). For simplicity, let the exterior 
medium be free space. 

From the surface-equivalence theorem of Schelkunoff [10], the 
interior field will vanish (without disturbing the exterior field) if 
we introduce the following surface-current densities 

(1 ) J = n x H 

. — s — 


( 2 ) 


M = E 
— s — 


x n 


on the closed surface S of the scatterer. (The unit vector n is directed 
outward on S.) In this situation, illustrated in Fig. lb, we may replace 
the scatterer with free space without disturbing the field anywhere. 

By definition, the incident field ( E . ,H . ) is generated by (J_. ,M . ) in 
free space, and the scattered field is: 11 


(3) 




(4) ^ > H - H, . 

When the surface current (J^Ms) radiates in free space, it generates 
the field (Es>Hs) in the exterior and (-Ei,-Hj) in the interior region. 
This result, illustrated in Fig. 1c, is deduced from Fig. lb and the 
superposition theorem. 

With the scatterer replaced by free space, we have noted in Fig. lb 
that the interior region has a null field. As shown in Fig. 2, we place 
an electric test source J_. in this region and find from the reciprocity 
theorem that 1 


(5) 


ft 


G* 


it - 2s • £t> ds + 


(J. • E t - M. • H t ) dv = 0 


where (^,Ht) is the free-space field of the test source. In words, 

Eq. (5) states that the interior test source has zero reaction with the 
other sources. This "zero-reaction theorem" was developed by Rumsey[6], 

Equation (5) is the integral equation for the scattering problem, 
and our objective is to use this equation to determine the surface- 
current distributions Jg and Kg. To accomplish this, we expand these 
functions in finite series so there will be a finite number N of unknown 


2 



Fig. la. The source ( J_-j ,M. ) generates the field (E_,H) 
with scatterer. 



Fig. lb. The interior field vanishes when the currents (J ,M ) 
are introduced on the surface of the scatterer. ~ s 


I \ 

K-e'.-h') \ 

\ ' \ 

, l ' 

( l 

FREE SPACE.' 

V ' 


Fig. lc. The exterior scattered field may be generated by 
(Jg.Mg) in free space. 
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expansion constants. Next we obtain N simultaneous linear equations to 
permit a solution for these constants. One such equation is obtained 
from Eq. (5) each time we set up a new test source. 





V>.\FREE SPACE 
x' 



FREE SPACE 


Fig. 2. An electric test source is positioned in 
the interior of the scattering region. 


The magnetic current M 5 vanishes if the scatterer is a perfect 
conductor. We assume a finite conductivity and use the impedance 
boundary condition: 


(6) = Z s ^ x S 


where Z s denotes the surface impedance. 

For 3-dimensional problems involving arbitrary scatterers, Js 
and Mg are functions only of the position on the surface of the scat- 
terer. If M. vanishes, Eqs. (5) and ( 6 ) yield 


(7) 


JJi 


4 


[E - (n 
L — m 


x HJ 
— m 


z s ] ds = 


ii 


E ds 
— m 


where (E ,H ) denotes the free-space field of test-source m. 

— m — m 

We represent the electric current distribution as follows: 

N 

(s) y i„i, 

n=l n ^ 

where the complex constants I are samples of the function J s< The 
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vector functions in are known as basis functions, subsectional bases, 

expansion functions or dipole modes. We employ expansion functions Jn 

and test sources J with unit current density at the terminals. 

— m 

From Eqs. (7) and (8) we obtain the simultaneous linear equations 
‘ 9 > j, W\i Wlthm=l,2, 3, -N 

where 

* r r 

0°) C mn = - J • [E - (n x HjZ c ] ds = - in * In ds 
mn „ — n L -m v -nr s m ~m — n 

(11) A m = ff J, • E ds = ff J • E. ds . 
m )). -l -m 

In Eqs. (10) and (11) the integrations extend over the region where the 
integrand is non-zero. For example, region n is that portion of the 
surface S covered by the expansion function J|>. Region m covers the 
interior test source Jxn- The reciprocity theorem relates the first 
and second integrals in Eq. (10). In the second integral, In is the 
free-space field generated by and the associated magnetic current 

For computational speed and storage, it will be advantageous to 
have a symmetric impedance matrix C m n. Furthermore, the test sources 
should be selected to yield a well-conditioned set of simultaneous linear 
equations. For these reasons, we employ test-sources Jjn of the same size, 
shape and functional form as the expansion functions In- Finally we 
position the interior test sources a small distance 6 from surface S 
and take the limiting form of the integrals as 6 tends to zero. 

The next section discusses the electric surface dipoles which are 
employed as test sources and expansion modes. 
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III. THE SINUSOIDAL SURFACE DIPOLES 


A planar surface dipole located on the yz plane is illustrated in 
Fig. 3a. This source is an electric surface-current density with height 
a and width b. The surface-current density is given by 


( 12 ) 


J = z cos 


^(z-z 2 ) \ 

l 2( z 2 _z i ) / 


for z 1 < z < z 2 


(13) 


( tt(z-Z 2 ) } 
l 2 l z 3' z 2 ) / 


for z 2 1 z 1 z 3 


As illustrated in Figs. 3b and 3c, the current density vanishes at the 
edges z = z i and z = z 3 , and is uniformly distributed in the transverse 
direction. The surface-current density and its slope are continuous 
across the terminal at z = z 2 « 

The sinusoidal surface dipole is a hypothetical source in free space. 
The current distribution on a rectangular plate is not sinusoidal. 

Figure 4 illustrates a surface V-dipole. Distance along the dipole 
arms is measured by the coordinates s and t with origin at the terminal 0. 
The surface-current density is 


(14) £ = -s cos | jj|— j on arm s 


(15) J. = t cos | on arm t. 

When the wedge angle ip is adjusted to 180°, the V-dipole in Fig. 4 
reduces to the planar dipole in Fig. 3a. 

Having defined the sinusoidal surface dipole, one is now in a 
position to explain its relevance. The dipole current distribution 
(Eqs. (14) and (15)) will be used as the basis functions ( Jji in Eq. (8)) 
for expanding the unknown current distribution induced on a conducting 
surface. Furthermore, surface dipoles will be employed as test sources 
with the reaction concept to solve the integral equation. 
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By superposition, the field of the surface dipole shown in Fig. 4 
is the sum of the field contributions from monopoles s and t. 


IV. THE IMPEDANCE MATRIX 

From the viewpoint of reaction, the complex number C mn in Eq. (10) 
represents the reaction between two sources. Let <m,n> denote the 
reaction between sources m and n, then Eq. (10) can be written as 


(16) 


C = -<m,n> = 

mn 


'm 


Ai 


L, ds 

— n 


It has been pointed out [6] that the reaction between two sources is 
related to the circuit parameters by 


(17) 


C = V I 
mn mn mm 


where V mn is the open circuit terminal voltage induced at m by source n, 

and I is the terminal current of source m when it transmits, 
mm 

Although the surface dipole described in Section III is a hypothetical 
source, it is useful to define its self- impedance with the induced-emf 
formulation: 


(18) 


mm 


mm 


rim 


C 

mm 



From Eq. (16), Eq. (18) yields 


(19) 


-1 


'mm 


mm 


m 


J 

— m 


E ds 
— m 


where Jm is the surface-current density of source m and Ejn is its free 
space electric field. The self-impedance of a center-fea, square 
surface dipole as a function of size is listed in Table I. 


The mutual impedance between two surface dipoles is defined by 


(20) 


mn 


= [ 

L I nn I 

mm nn J 


'm 


4 


E ds 
— n 
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TABLE I 


Self Impedance of Center-Fed Planar 
Surface- Dipole Shown in Fig. 5 


a/ A 

R n 

x n 

0.2 

11.48 

-69.76 

0.3 

23.98 

-35.26 

0.4 

38.80 

-15.36 

0.5 

52.98 

- 4.52 



Fig. 5. Surface Dipole with J = z cos (irz/a) 
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Figure 6 illustrates a pair of center-fed, coplanar surface dipoles, and 
Table II lists the mutual impedance Z, 2 . Here s and s specify the 
relative position of the dipoles. z y 


V. THE EXCITATION COLUMN 

The complex quantities Am in Eq. (11) form the excitation column 
in the matrix equation C mn I n = Am. Physically, A m is the reaction 
between the impressed source and test dipole m. From Eq. (11), 

A^ is given by 


( 21 ) 


m 


JJ 


m 


in ' E, ds 
— m —i 


The above integral requires numerical integration over the source m. 

If the source Jj 1S located at a great distance from test dipole m, the 
incident field (ji • ,H^ ) may be regarded as a plane wave with 

jk(x'sin9'. cos<f>'. + y'sinel sin<}>'. + z'cose'.) 

(22) E. = ^ e 11 11 1 


where (r 1 , e-} , 4 >-} ) are the spherical coordinates of the source and Eq is 
the incident electric field intensity at the coordinate origin O'. 
Figure 7 illustrates an incident plane wave illuminating an electric 
surface dipole located on the y'z' pjane with height a and width b. 

For the surface-current density J. “ z‘ cos(irz , /a)» Eq. (21) can 
readily be evaluated to yield 


(23) 


m 



z' ) 27rab 


sin(X i ) cos(Y.j/2) 

X, (v? - u 2 ) 


where 

X. = 0.5 kb sin el sin$l 
Y i = ka cosej. 
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TABLE II 

Mutual Impedance of Center- Fed Coplanar Surface Dipoles Shown in Fig. 6 

a/X = 0.5 b/x = 0.25 


0.75 

1 .076- j 7.913 

-2 . 925- j 5.950 

-5 . 968+ j 1.761 

0.8871+j 7.019 

0.50 

24.43 +j 5.997 

10.04 -j 9.321 

-9 . 568- j 6.65 

-7.194 +j 8.971 

0.25 

53.10 +J55.87 

29.54 -j 9.070 

-9 . 468- j 19.33 

-15.74 +j 6.426 

0.0 

67.09 + j 13.23 

39.07 - j 22 .49 

-8.659-j26.75 

-19.32 +j 4.263 


0.0 

0.25 

0.50 

0.75 



l_l 

*-b-> 


Fig. 6. Coupled surface dipoles 
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X 


Fig. 7. A plane wave (E.,H.) illuminates an electric 
surface dipoleT 1 
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VI. FAR-FIELD RADIATION AND SCATTERING 


The field scattered by a perfectly conducting body may be generated 
by the electric surface-current density ^ in free space. To obtain the 
total field one adds the free-space field (Ej ,Hj ) of the electric source 

-i- 


Consider an electric surface dipole with current density 
Js = z' cos(-rrz7a) located on the y'z 1 plane as shown in Fig. 8. From 
reciprocity, the free-space electric field generated by this source at 
a distant point (r‘ ,e;,<f>;) is 


(24) 


E s = e' 


• sin(X_) 

OifiiL a b §_ 

s 2 aD X„ 


cos(Y $ /2) 

1 ^ 7 ) 


-jkr' 


sine; 


where 


X s = 0.5 kb sine; sin<j>; 

Y = ka cose' . 
s s 

a /v 

The, e- and <f>- components of the scattered field with respect to the 
reference coordinate system 0 can be obtained easily via an appropriate 
coordinate transformation. 

In plane-wave scattering problems, one is often interested in the 
echo area a defined as follows 

|E S | 2 

(25) a = lim 47rr 2 

r-x» [E 1 | 2 


where E 1 is the incident electric field intensity. 

In antenna and radiation problems, one is interested in the power 

gain: 


4tt r 2 |E.| 2 

(26) Gain = 5 

n | V | 2 G 


where V is the terminal voltage of the antenna, G is the terminal con- 
ductance of the antenna, and n is the intrinsic impedance of free space. 


14 ' 


z 



/ 


Fig. 8. A surface dipole radiates in free space. 


VII. NUMERICAL RESULTS 

Figure 9 presents the backscattering echo area of a square plate with 
perfect conductivity for the broadside aspect.. In the reaction cal- 
culation, the plate is divided into cells, and overlapping current 
modes were employed as illustrated in Fig. 10. In this case the trans- 
verse current is neglected and '45 modes were used for the current 
distribution. Useful results can be obtained with as few as one mode 
per square wavelength of surface area. For comparison. Fig. 9 also 
shows the experimental measurements of Kouyoumjian [11]. 

The magnitude and phase of the current density induced on a 
perfectly-conducting rectangular plate are illustrated in Figs. 11 and 
12. Figures 13 through 16 show the normalized backscatter cross-section 
of a rectangular plate. Figures 17 and 18 show the normalized backscatter 
cross-section of a corner reflector. The title of each figure gives the 
echo area at the broadside aspect in terms of dB = 10 log (a/ A 2 ). 
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Fig. 9 '. Backscatter cross-section of perfectly conducting 
square plate for the broadside aspect. 


h 

h 

ii 



-II- 

-12- 

-13- 

-14- 

-15- 

-6- 

— 7 — 

— 8 — 

— 9 — 

-10- 

— 1 — 

-2- 

-3- 

-4- 

-5- 


Fig. 10. Electromagnetic modeling of plate. 
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20 , 



- 0.5 - 0.25 0 0.25 0.5 

y/x 


Fig. 11. Amplitude of surface-current density induced on a 
perfectly-conducting rectangular plate for a 
plane wave incident at broadside. 



y/x 


Fig. 12. Phase of surface-current density induced on a 
rectangular plate for a plane wave incident at 
broadside. 
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REACTION 


MEASURED 


H = 3X W*2X 



9 (DEGREES) 


Normalized backscatter cross-section in the yz plane 
of a rectangular plate, a (e,<|>) = 15.25 dB at (90°, 90°) 


REACTION 

MEASURED 

H = 2X W=3X 





90 

9 ( DEGREES ) 


Fig. 14. Normalized backscatter cross-section in the yz plane 

of a rectangular plate. a Q .,(e,<j>) = 15.23 dB at (90°, 90°) 







Fiq. 15. Normalized backscatter cross-section i'n the xy plane of 
a rectangular plate. cj QQ (e,<j>) = 15.25 dB at (90°, 90°). 



ig. 16. Normalized backscatter cross-section in the xy plane of 
a rectangular plate. = 15.23 dB at (90°, 90°). 


REACTION 

: MEASUREO 

H = X W = 0.5X a = 1 30 


(DEGREES) 


Normalized backscatter cross-section in the xy plane of 
a corner reflector. a QQ (e,<|>) = -0.38 dB at (90°, 90°). 

0 O v 


REACTION 

MEASUREO 

H = X W=0.5X a = 130 


T 

O <r 98< <1B1 


8 (DEGREES ) 


Normalized backscatter cross-section in the yz plane of 
a corner reflector. = - .38 dB at (90°, 90°). 


m 



Figures 20 through 23 show the Emplane gain of the corner-reflector 
antenna illustrated in Fig. 19. Figures 24 through 27 show the H-plane 
gain of the same antenna. For comparison. Figs. 20 through 27 include 
experimental measurements obtained by Melvin Gilreath at NASA Langley 
Research Center. In the experimental measurements the receiving antenna 
was linearly polarized in the theta direction. Similarly, the cal- 
culated gain is based on E . The dipole length is x/2 and the radius 
is 0.005X. 6 


In the reaction calculation, only vertical modes were employed to 
approximate the current distribution. The number of modes used to obtain 
the results given in Figs. 11 through 27 are listed below. In each case, 
the matrix size is equal to the number of modes. 


VIII. 


Figs. 

Number of 

11,12 

45 

13,14 

55 

15,16 

75 

17,18 

30 

20-27 

61 

CONCLUSIONS 



The reaction concept and Galerkin's method are used to develop a 
new formulation for perfectly-conducting antennas and scatterers. Numer- 
ical results are presented for scattering and radiation from rectangular 
plates and corner reflectors. The results show general agreement with 
measurements. 


The techniques can be applied to surfaces with finite conductivity 
and arbitrary shape. For arbitrary polarization and aspect, 24 modes, 
per square wavelength are adequate. This is a significant improvement 
over previous techniques. 
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-10 -20 


-20 


-10 
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Fig. 20. Relative gain in the E-plane of a corner-reflector antenna. 
G(e,<j>) = 4.31 dB at (90°, 90°). 
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Fig. 21. Relative gain the the E-plane of a corner-reflector antenna. 
G ( 0 , 4 > ) = 4.05 dB at (90°, 90°). 
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Fig. 22. Relative gain the the E-plane of a corner- reflector antenna. 
G(e,<(>) = 7.48 dB at (90°, 90°). 
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Fig. 23. Relative gain in the E-plane of a corner-reflector antenna. 
G(e,<|>) = -1.06 dB at (90°, 90°). 
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• ••REACTION S = 0.75 X a = 90° 



' Fig. 24. Relative gain in the H-plane of a corner-reflector antenna. 
G(e,<j>) = 4.31 dB at (90°, 90°). 
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Fig. 26. Relative gain in the H-plane of a corner-reflector antenna. . 
G(e,<j>) = 7.48 dB at (90°, 90°).. 
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I 


9 © 9 REACTION 
MEASURED 



Fig. 27. Relative gain in the H-plane of a corner-reflector antenna. 
G(e,<j>) = -1.06 dB at (90°, 90°). 
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